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第 一 性 原理 在 热 中 子 散 射 截 面 计 算 中 的 应 用 
Eu, 张 信 一 ŽE, MER, Eg) 
Cl. 西北 核 技 术 研 究 所 ， 西 安 ，710024) 
摘 要 : 第 一 性 原理 被 应 用 在 核 工程 领域 关注 的 热 中 子 散射 材料 的 热 中 子 截面 计算 中 , 本 文 以 Al 和 Bi 金属 的 
热 中 子 截面 制作 为 例 ， 分 别 基 于 第 一 性 原理 冷冻 声 子 法 和 密度 泛 函 微 扰 法 ， 采 用 VASP 结合 PHONONY 程序 
计算 了 它们 的 声 子 色 散 关 系 和 声 子 态 密度 ， 基 于 NJOY 程序 ， 在 LEAPR 中 添加 Bi 相干 散射 处 理 ， 制 作 了 Al 
和 Bi 的 热 中 子 散射 截面 库 ， 结 果 表 明 : 对 于 Al， 相 较 于 冷冻 声 子 法 ， 密 度 泛 函 微 扰 理论 的 声 子 态 密度 得 到 的 
热 中 子 散射 截面 与 ENDF8.0 符合 更 好 ， 对 于 Bi， 采 用 密度 泛 函 微 扰 法 消除 了 冷冻 声 子 法 的 虚 频 现象 ， 得 到 的 
热 中 子 散 射 结果 与 实验 值 符合 较 好 。 本 文 从 探索 材料 内 部 特征 机 理 的 角度 提出 了 一 种 更 基础 、 更 有 预见 性 的 
制作 热 中 子 截面 的 方法 ， 为 研究 新 型 反应 堆 的 核 材 料 热 化 机 理 莫 定 了 理论 基础 。 
关键 词 : 第 一 性 原理 ， 声 子 态 密度 ， 热 中 子 散 射 截面 ， 核 数据 制作 
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反应 堆 工 程 领域 对 计算 机 模拟 工具 有 着 广泛 的 应 用 , 近年 来 ， 国外 发 展 的 先进 模拟 与 仿真 
技术 依托 高 性 能 计算 能 力 ， 以 “第 一 性 原理 ”的 物理 学 为 基础 建立 模型 ， 取 代 经 验 公 式 ， 在 虚拟 
世界 中 重建 现实 世界 的 实体 与 系统 ， 实 现 对 复杂 系统 的 预测 性 模拟 ， 已 在 材料 研发 、 气 候 预 测 
等 领域 取得 了 巨大 进展 与 效益 。 核 科学 技术 研究 和 核 工 程 设计 中 ， 为 了 提高 核 设 计 的 精确 度 ， 
可 从 两 方面 入 手 : 一 方面 是 努力 改进 核 设计 的 计算 模型 和 计算 方法 ， 以 提高 计算 的 精确 度 ， 另 
一 方面 是 提高 核 数据 的 精确 性 。 对 于 核 工程 技术 人 员 来 讲 ， 正 确 地 了 解 和 使 用 核 数据 是 非常 重 
要 的 ， 因 为 它 是 获得 正确 计算 结果 的 前 提 和 基础 。 

热 中 子 反 应 堆 设计 计算 使 用 的 核 材料 截面 , 尤其 是 慢 化 剂 材 料 的 热 中 子 散射 数据 对 热 中 子 
反应 堆 的 核 设计 和 中 子 学 分 析 、 热 中 子 过 滤器 、 冷 中 子 源 设 计 等 工作 至 关 重 要 。 热 能 区 中 子 的 
散射 截面 不 仅 与 中 子 的 能 量 相 关 ， 还 与 发 生 散射 介质 的 温度 、 物 理 和 化 学 性 质 相 关 ， 强 烈 依 赖 
材料 本 身 的 声 子 态 密度 的 性 质 , 以 往 的 常用 慢 化 剂 材料 的 声 子 态 密度 大 多 数 基于 上 世纪 六 十 年 
代 的 简化 声 子 模型 叫 ， 在 引入 较 多 近似 的 条 件 下 才能 得 到 热 中 子 散射 数据 ， 给 反应 堆 临 界 和 中 
子 能 谱 计 算 引入 了 较 大 误差 , 近 十 年 来 国内 外 陆续 基于 第 一 性 原理 和 分 子 动力 学 的 理论 进行 了 
多 种 材料 的 热 中 子 散 射 数据 产生 与 计算 方法 的 研究 ， 并 持续 更 新 国际 核 评价 数据 库 C53。 早 期 
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的 ENDF/B-VL8 只 有 十 五 种 热 中 子 散射 材料 ，ENDF/B-VIL0 有 二 十 种 热 中 子 散 射 材料 ， 
ENDF/B-VIL.1 扩充 到 二 十 一 种 热 中 子 散射 材料 ， 最 新 的 ENDF/B-VILO 更 是 提供 了 二 十 四 种 
热 化 材料 的 三 十 四 项 热 散 射 数据 ， 同 一 材料 不 同 结构 的 情况 得 到 了 有 效 处 理 ， 比 如 石墨 ， 数 据 


库 就 给 出 了 三 种 不 同 孔 际 率 的 热 散 射 数据 , 扩充 的 这 些 热 中 子 散射 数据 库 大 部 分 是 基于 第 一 性 
原理 和 分 子 动力 学 模拟 得 到 的 声 子 态 密 度 完 成 的 。 热 中 子 材 料 声 子 态 密度 的 计算 需要 了 解 原子 
系统 的 动力 学 ， 这 包括 唱 格 结构 和 由 声 子 色散 关系 和 偏振 矢量 表示 的 振动 类 型 ， 


这 种 振动 可 以 


以 声 子 态 密度 的 形式 表示 ， 而 传统 的 声 子 态 密度 是 通过 建立 原子 间 相互 作用 的 解析 模型 、 评 估 


原子 间 的 力 常 数 、 构 造 动力 学 矩阵 和 对 角 化 来 计算 的 ,动力 矩阵 计算 中 的 力 常 数 可 以 用 不 同 的 


方法 来 估计 ， 一 般 它 们 来 自 热力 学 性 质 ， 比 如 比热容 或 者 压缩 率 数据 ， 此 外 ， 还 可 以 利用 非 弹 
性 散射 技术 从 声 子 色散 关系 的 实验 测量 中 推断 出 声 子 色散 关系 , 这 两 种 方法 基本 上 都 是 半 解 析 
关系 式 和 实验 数据 的 拟 合 ， 有 两 个 主要 缺陷 ， 第 一 ， 它 不 是 可 预测 的 ,产生 的 原子 力 和 常数 的 色 
散 关系 是 由 实验 数据 推测 的 ， 第 二 ， 结 果 并 不 是 唯一 的 ， 有 可 能 被 其 他 动力 学 模型 奉 代 。 相 比 
于 半 经 验 关 系 式 产生 声 子 态 密度 方法 ， 第 一 性 原理 方法 更 贴近 物理 本 质 ， 无 需 拟 合 实验 数据 ， 


它 将 多 元 素 的 原子 体系 当 作 由 原子 核 和 电子 构成 的 多 粒子 体系 ,针对 不 同 的 唱 格 结构 ， 基 于 量 
子 力学 的 基本 原理 ,在 不 引入 任何 实验 参数 的 条 件 下 对 多 原子 体系 进行 处 理 ， 基 于 严格 和 准确 


问题 的 求解 归结 为 


的 密度 泛 函 理 论 ， 将 固体 抽象 为 具有 平移 周期 性 的 理想 晶体 ， 认 为 原子 围绕 其 平衡 位 置 振荡 ， 


BE 子 在 周期 性 势 场 中 的 运动 的 求解 。 利用 第 一 性 原理 计算 材料 声 子 态 密度 ， 


较 之 前 凭借 经 验 公式 和 简单 的 动力 学 模型 的 计算 ， 具有 明显 的 通用 性 ， 它 允许 了 更 复杂 、 更 系 


统 的 多 原子 动力 学 模型 ， 为 热 中 子 散 射 截面 的 计算 莫 定 了 重要 的 理论 基础 。 


为 满足 热 中 子 反应 堆 对 热 中 子 散射 数据 的 迫切 需求 , 针对 目前 热 中 子 散 射 数据 产生 与 处 理 


技术 面临 的 问题 ， 本 文 开 展 了 第 一 性 原理 的 热 中 子 散射 截面 计算 研究 。 本 文 


要 分 为 以 下 几 个 


部 分 : 第 一 部 分 简要 概述 了 热 中 子 散射 理论 ; 第 二 部 分 介绍 了 第 一 性 原理 Al 和 Bi 金属 热 中 子 散 


射 数 据 的 计算 方法 和 建 模 ， 第 三 部 分 给 出 Al 和 Bi 金属 热 中 子 散 射 数据 的 计算 结果 并 讨论 分 析 ; 


第 四 部 分 给 出 了 本 文 的 结论 。 
1 热 中 子 散 射 理 论 
在 热 中 子 反应 堆 中 , 对 于 能 量 低 于 1 eV 的 中 子 , 由 于 中 子 能 量 与 散射 核 上 


的 热能 是 可 比较 的 ， 


不 能 简单 认为 靶 核 是 静止 的 ， 这 种 情况 下 ， 中 子 很 容易 与 物理 离子 发 生 非 弹 


性 


nbd, ENT 


Fa AR Sree WAR eae REEE, 通过 热 中 子 非 弹性 散射 截面 表示 ,同时 由 于 较 低能 


量 中 子 的 德 布 罗 意 波长 与 物质 分 子 或 晶体 内 书 


Xr 


的 间距 相差 不 大 , 与 不 同 核发 生 散射 的 中 子 间 可 
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能 发 生 王 涉 ， 即 布拉格 效应 ， 通 过 相干 弹性 散射 截面 来 表示 四。 热能 区 的 中 子 散射 截面 不 单纯 
是 与 中 子 能 量 有 关 ， 还 与 散射 介质 的 温度 及 物理 、 化 学 性 质 有 关 ， 它 反映 的 是 材料 自身 的 声 子 
态 密度 特征 ， 因 此 必须 要 研究 散射 物质 的 声 子 振动 特性 。 

对 非 弹性 散射 截面 , 包括 相干 和 非 相 干 非 弹性 散射 , 对 所 有 物质 都 重要 , 用 热 散 射 律 表示 。 
非 弹性 散射 截面 一 般 由 其 双 微 分 散射 截面 推导 而 来 中 : 


2 
o;,(E,T) = [| <A OME (1) 


其 中 E 和 E' 是 实验 系 中 入 射 和 散射 中 子 的 能 量 ，Q 是 实验 系 下 的 散射 角度 ，T 是 温度 ， 单 
位 是 eVv， 为 了 计算 得 到 非 弹 性 散射 截面 ， 一 般 利 用 NJOY 程 序 的 LEAPR 模 块 和 THERMR 模 块 ， 
NJOY 程 序 的 计算 主要 基于 “不 相干 近似 ”BlJ， 认 为 : 


2 ' 
: —_ a Ë Iont Fne|8,(0t, B) 


O con ERRA AIA PON BC Ope 是 材料 的 束缚 核 的 非 相 干 散射 截面 ，a 为 动量 转移 


量 ，B 为 能 量 转 移 量 ， 表 达 式 分 别 如 下 : 


jE+E-2YE'E 0080 _ ie 


AkT ~ 2MkT 
Ss:(0,B) 为 不 考虑 相互 作用 的 自 散射 律 ， 可 以 写成 
S$,(0,P)=7| ee’ Od © 
KP, 7 QUE AD), NAT) K y 表示 为 
Hi) =af* LE j-e © 


»2Bsinh (8/2) 


其 中 的 p(B) EU BASAL PAS, 并且 是 归 一 化 的 ， 也 是 NJOY 程序 唯一 需要 输入 的 
外 部 量 ， 也 是 本 文采 用 第 一 性 原理 计算 的 重要 输出 值 ， 声 子 态 密度 的 结果 直接 影响 热 中 子 散 射 
数据 的 质量 和 准确 度 。 


当 Q AB ABEAR (a> @,,, 


BI> B,, )， 即 在 较 高 的 入 射 能 量 的 时 候 ，Q AB 已 经 超出 


[Some 


了 S(@ 户 范围 , RORY at ie BES] A RIN Te AEA (SCT), FAH T RER ERGIT KAA F At asl: 


_ (appr IA 


= 4aT G(T) t3] 
Ss (æ, B, T) = == 7 
aa FR o 
4ra = 
T 
BEMIS EI AA EERE PRA R A 
EVE (a-|B)’ T B+ 
oT E= 7 T r 
l Ag 7 ae z] 4da To 2 ®) 
其 中 有 效 温 度 是 由 声 子 谱 计 算得 到 的 : 
1 lre ho 
Ten = 5 J hocoh 22 poao (9) 


对 于 由 相干 散射 体 组 成 的 物质 中 ， 零 声 子 项 给 出 了 相干 散射 函数 ， 对 应 的 相干 散射 截面 为 


Onn, (E, u) = A > fer" b(u > L) (10) 


E,<E 


RP, 是 跟 晶 格 结构 因子 有 关 的 一 个 系数 : 


2mh? 
S E of (11) 
五 (7) 是 每 个 原子 的 相位 ， 表 达 式 如 下 : 
2 
Fo 了 Fei (12) 
对 (10) 式 进行 角度 和 能 量 的 积分 ， 得 
coh 一 > n (13) 
E<E 
2,2 
已 = (14) 


其 中 ， 丈 是 倒 格 子 壳 矢量 长 度 ，m 是 中 子 质量 


不 同 唱 格 结构 唱 体 的 结构 因子 计算 是 不 同 的 ，NJOY 的 LEAPR 模 块 中 仅 内 置 了 面 心 立方 


FCC (Face-Centred Cubic)、 体 心 立 方 BCC (Body-Centred Cubic ) 和 六 角形 HCP (Hexagonal Close 


Packed) 的 结构 ， 对 


FH HSA BS 
LEAPR 模 块 ， 比 如 对 于 Bi 金属 的 相干 弹性 散射 ， Bi 金属 属于 


结构 不 符 的 结构 就 需要 用 户 自己 输入 结构 因子 ， 并 重新 编译 
BA mA, 单元 晶 胞 中 含有 两 个 原 


子 ， 菱 唱 角 度 为 57.23*， 在 298.0 K 温 度 下 Bi 原子 的 位 置 参 数 z=0.2271， 倒 格子 壳 矢 量 可 以 表示 
为 
t=1b +1,b,+Lb, (15) 
其 中 4，1， ，4 是 晶 面 指数 ， 为 表征 晶 胞 位 置 的 参数 ，b ”1, 为“， 万 “3 是 坐标 矢量 ， 表 达 
式 如 下 
hap, ee: E a 
a, ʻa, XA; ay Ba, Cy 
Ja E a (16) 
; aʻa, Xa; Ba, Cy 
a, Xa, 1 | . thx 
b, = 27 =27( yt+—z) 
a A, XQ, V3ay H 


Gyy，Cp 为 对 应 六 角形 坐标 系 
构 包 含 两 个 Bi 原子 ， 
因子 


Fol=je 
2 计算 方法 与 建 模 


Fete 


Bi 的 晶 胞 参数 ，a， =4.557A, ，c，= 11.862 A ， 原 胞 结 


i272(1+1,+1,) 


分 别 位 于 坐标 (Z,Z,Z) 和 (-Z,-Z,-Z) 处 。 由 公式 (12)， 可 以 得 到 Bi 的 晶 格 结构 


q 2h th +h) 2 (17) 


=4cos [2zz(0 +1, +1,)| 


本 文 首 先 研 究 了 第 


© 根据 畸变 品格 和 理想 晶 格 的 能 量 差 来 计算 出 位 移 振幅 


计算 出 声 


性 原理 声 子 态 密度 的 产生 方法 , 可 以 分 为 两 类 , 第 一 


类 为 直接 计算 法 ， 


子 态 密度 ， 称 为 冷冻 声 子 法 ; 第 二 


函数 ， 进 而 通过 移 位 原子 得 到 力学 矩阵 ， 
方法 为 间接 方法 ， 通 过 描述 价 电子 密度 对 周期 品 


格 的 扰动 响应 得 到 声 子 态 密度 结果 ， 


称 为 密度 泛 函 微 扰 理论 法 。 和 常用 的 软件 有 CASTEP 程 序 、 


VASP 程 序 、PHONON 程 序 和 PHONONPY 程 序 ， 


其 中 VASP 和 商业 程序 PHONON 的 求解 是 采用 


的 是 基于 密度 泛 函 微 扰 理论 


的 冷冻 声 子 FP(Frozen Phonon, FP) 直 接 法 ， 而 CASTEP 采 月 


简称 DFPT(Density Function Perturbation Therory) 方 法 ， 


PHONONPY 通 


过 脚本 的 方式 既 可 以 采用 FP 法 又 可 以 基于 DFPT 方 法 。 


的 方法 ， 


本 文采 用 VASP 程 序 结合 开源 程序 


如 图 1 所 示 为 第 一 性 原理 冷冻 声 子 FP 法 计生 


H 


ETSK 


度 的 流程 图 , 首选 选择 要 计算 的 唱 胞 ， 


产生 它 的 超 品 胞 组 成 结构 ,通过 
过 力 常 数 矩 阵 的 傅 里 叶 变 换 得 到 动力 学 矩阵 ， 


动 


再 通过 


寺 对 称 性 移 位 原子 , 计算 所 有 原子 的 受 力 情况 , 构建 力 常数 矩阵 ， 


力学 矩阵 对 角 化 计算 材料 的 声 子 色散 


曲线 ， 分 析 材 料 的 声 子 色散 曲线 ， 得 到 材料 的 声 子 态 密度 ， 利 用 得 到 的 声 子 态 密度 就 可 以 得 到 
材料 的 热 中 子 散射 数据 和 相关 的 热力 学 性 能 参数 ， 它 是 一 种 直接 的 方法 ， 基 于 原子 从 平衡 位 置 
位 移 的 总 能 量变 化 或 力 计算 的 方法 ， 只 要 原子 的 受 力 可 以 导出 ， 就 可 以 根据 力学 常数 矩阵 得 到 
相应 的 声 子 态 密 度 ， 但 是 该 方法 的 缺点 是 必须 要 考虑 超级 唱 胞 ,这 对 于 单元 与 单元 之 间 扩 展 较 
小 的 简单 晶体 结构 是 不 太 适 用 的 。 而 密度 泛 函 微 扰 理 论 DFPT 法 为 间接 计算 法 ， 它 是 对 平衡 品 
格 几 何 进 行 扰动 后 ， 通 过 微 扰 理 论 分 析 计 算得 到 力 常数 的 解析 方法 ，DFPT 不 同 于 冷冻 声 子 方 
法 ， 其 主要 的 原理 是 通过 计算 晶体 对 外 界 的 能 量变 化 的 响应 从 而 计算 晶体 的 动力 学 性 质 ， 该 方 
法 不 用 扩 胞 ， 也 不 需要 波 矢 正 交 与 晶 胞 边界 ，DFPT 为 复杂 晶体 的 唱 格 动力 学 分 析 提 供 了 有 力 
的 计算 手段 ， 被 认为 是 相当 精确 、 有 效 和 可 靠 的 方法 ， 然 而 它 需 要 高 度 专业 化 的 第 一 性 原理 计 
算 程 序 ， 并 且 还 需要 基于 原始 程序 对 其 进行 相关 的 改动 和 补充 ， 这 也 是 本 文采 用 开源 程序 
PHONONPY 的 初 囊 。 本 文 将 分 别 对 比 这 两 种 方法 在 晶体 Al 和 Bi 声 子 态 密度 计算 中 的 差异 ， 蝇 


体 Al 和 Bi 的 唱 格 结构 如 图 2 所 示 ， 其 中 Al 属于 典型 的 面 心 立 方 结构 ，Bi 属 于 葵 方 晶 系 。 
选择 晶 胞 结构 
一 
产生 超 晶 胞 
+ 


| 通过 对 称 性 移 位 原子 ， 计 算 所 有 原子 的 受 力 情况 


$ 


| 构建 力 常数 矩阵 | 


ý 


| 通过 力 常 数 矩 阵 的 傅 里 时 转换 计算 动力 学 矩阵 | 


通过 动力 学 矩阵 对 角 化 计算 声 子 色 散曲 线 


| 计算 声 子 态 密度 | 
y 


| 计算 热 中 子 散射 截面 和 热力 学 性 质 等 


图 1 第 一 性 原理 冷冻 声 子 FP 法 计算 声 子 态 密 度 流程 图 
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(a) Al 原子 (b) Bi 原子 
图 2 Al 和 Bi 的 晶 胞 结构 


在 NJOY 制 作 热 散 射 律 的 过 程 中 需要 注意 一 些 操作 细节 ， 如 图 3 所 示 为 NJOY 制 作 ACE 各 式 
的 热 中 子 散射 律 数据 库 的 具体 步骤 , 其 中 MODER 模 块 是 将 ENDF/B 库 中 原始 的 核 素 数据 库 进 行 
十 进 制 到 二 进 制 的 转换 ,因为 计算 机 便于 运行 二 进 制 文件 , 然后 RECONR 模 块 对 原始 核 素 库 的 
截面 数据 进行 共振 重 造 ， 重 建 误差 设 为 0.001， 紧 接着 BROADR 模 块 对 其 进行 多 普 勒 展 宽 ， 注 
化 误差 设 为 0.001， 然 后 将 加 工 后 的 快 中 子 数据 提供 给 THERMR 模 块 ， 与 此 同时 ，LEAPR 模 块 
将 加 工 好 的 ENDFB 格 式 的 热 中 子 散射 律 数据 给 THERMR 模 块 ，THERMR 模 块 负责 热 区 散射 矩 
阵 等 的 计算 ， 然 后 将 生成 的 点 截面 库 PENDF 格 式 的 数据 库 送 给 ACER 模 块 ， 制 作成 MCNP 软 件 
的 ACE 格式 的 数据 库 。ACER 模 块 最 后 生成 两 个 数据 文件 : 一 个 是 新 制 的 MCNP 中 子 截面 库 文 
件 , 放 到 MCNP 截 面 库 文件 夹 内 ; 另 一 个 是 路 径 文件 , 用 来 添加 到 MCNP 软 件 的 索引 文件 XSDIR 
里 面 , 以 供 MCNP 程 序 计算 时 调用 相应 截面 库 。 整个 流程 中 最 重要 的 是 LEAPR 模 块 和 THERMR 
模块 ,多 普 勒 展 宽 等 对 其 影响 很 小 ， 主 要 是 为 了 提供 一 些 基本 的 散射 截面 等 数据 ， 而 热 区 中 子 
在 核 素 内 的 散射 情况 由 LEAPR 模 块 的 散射 律 数 据 决 定 。 


BROADR 


LEAPR /TSL 库 (ENDF/B 格式 ) 


i 
F 


图 3 NJOY 制作 ACE 格式 热 散 射 截面 流程 图 


目前 NJOY 在 制作 热 散 射 律 过 程 中 ,常见 慢 化 剂 材料 的 MAT 号 和 相关 的 反应 道 MT 参 照 早期 
ENDF/B VI 格式 的 数据 库 定义 的 材料 及 热 截 面 号 码 ， 人 然而， 早期 ENDF/B VI 格式 的 数据 库 对 慢 
化 剂 材 料 的 定义 不 是 很 完善 ， 只 对 常用 慢 化 剂 材料 进行 了 研究 ， 导 致 NJOY 对 新 型 慢 化 材料 的 
识别 (MAT 号 码 ) 还 有 限制 ， 本 文 参照 ENDF/B VI 现 有 格式 的 规则 ， 针 对 Al 和 Bi， 重 新 定义 热 


pp 


[Some 


的 MAT 号 和 热 截 面 MT 号 码 ， 以 供 NJOY 的 LEAPR 和 THE 


中 子 散 射 得 


除了 材料 本 身 的 声 子 态 密度 ，NJOY 程 序 还 需要 输入 材料 的 束缚 核 截 面 等 基本 参数 ， 


RMR 使 用 。 


表 1 


列 出 了 Al 和 Bi 束缚 核 的 A 及 不 同 反 应 截面 的 值 ， 可 以 看 出 ， 对 于 Al 和 Bi 这 种 金属 晶体 ， 相 干 弹 
性 散射 是 比较 严重 的 ，Al 是 FCC 结 构 ，NJOY 程 序 中 内 置 了 相干 散射 计算 功能 ， 但 Bi 计算 时 需 
要 重新 构 置 晶 格 学 结构 因子 ， 并 且 添 加 到 THERMR 模 块 中 。 
表 1 Al 和 了 Bi 核 的 不 同 反应 的 束缚 截面 
原子 A O „p /barn O „on /barn Onc /barn O, /barn O fee barn 
7A] 26.75 0.231 1.495 0.0082 1.503 1.3966 
20?9Bi 207.2 0.0338 9.148 0.0084 9.130 9.0683 


3 计算 结果 与 讨论 
3.1 品格 参数 对 比 
本 文采 月 


的 厦 势 ， 利 用 VASP 程 序 包 


昌 VASPL0 进 行 AL 和 Bi 的 几何 结构 优化 ， 选 取 Pearson 数 据 库 中 的 实验 值 作为 几何 优 
化 的 初始 参数 ， 模 拟 唱 格 参数 和 原子 位 置 ， 分 别 选 取 广 义 梯度 近似 GGA 和 
的 结构 优化 选项 将 结构 松弛 到 它 的 最 低能 


局 域 密度 泛 函 LDA 


i=) 
里 


状态 ， 初 始 参 数 选取 


Pearson 数 据 库 中 的 实验 值 ,对 品格 参数 和 原子 位 置 进行 模拟 , 电子 能 


ee AC Np HEI A 1x10 eV, 


进行 
平面 波 截断 能 为 400 eV， 计 算 结 果 见 表 2 所 示 ， 从 结果 中 可 以 看 出 ， 


与 局 域 密度 泛 函 LDA 相 比 ， 


广义 梯度 近似 GGA 方 法 优化 得 到 的 唱 格 参数 更 接近 Pearson 数 据 库 
关 Al1 和 Bi 声 子 态 密 度 的 计算 就 采用 GGA 质 势 优 化 得 到 的 结构 进行 


表 2 Al 和 Bi 晶 格 参数 的 对 比 
GGA 
误差 
0.2% 
0.4% 
0.4% 


材料 方法 


值 
4.05 
4.57 
11.90 


Al a=b=c 


i a 
Bi 
c 


3.2 声 子 色散 关系 和 声 子 态 密 度 对 比 
日 优化 好 的 Al 和 Bi 的 晶 格 结构 ， 本 文 利用 了 灵活 性 更 强 的 
计算 了 Al 和 Bi 的 声 子 态 密度 , 交换 关联 函数 使 月 


AA 


VASP 3K (4 
eV, it SORA ATA 
网 格 Monkhost-Pack 格 子 中 进行 ， 本文 分 别 采 


利用 冷冻 声 子 FP 法 和 密度 泛 函 微 扰 理 


HGGA-PBE, 截断 
E 超 晶 格 方法 ， 超 晶 胞 大 小 为 4x4x4，Al 和 Bi 的 布 量 
] 了 FP 法 和 DFPT 法 计算 了 Al 和 Bi 


E 论 DFPT 法 得 到 Al 声 子 色 散 关 系 在 布 里 洲 


1 的 实验 值 ， 接 下来， 本文 


1T 0 


LDA 


Pearson 


-1.2% 
1.1% 
2.9% 


4.04 
4.56 
11.86 


PHONONPY 开 源 软件 结合 
最 终 选 取 为 400 
EE 渊 区 积分 在 3x3x3 的 KK 
的 声 子 色散 关系 ， 
区 不 同方 向 的 计 


ob E 


用 里 


202306.00654v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


算 结果 如 图 4、 图 5 所 示 ， 较 低 的 分 支 为 声学 项 ， 较 高 的 分 支 为 光学 项 。 从 图 中 可 以 看 出 ， 对 于 
Al 晶 体 而 言 ，DFPT 和 FP 结果 基本 类 似 ， 均 可 以 用 来 产生 合理 的 声 子 态 密 度 ， 进 而 得 到 热 中 子 
散射 数据 ， 对 Bi 晶体 ， 用 FP 方法 得 到 的 色散 关系 在 伽 马 点 周围 包含 一 些 虚 频 ， 这 可 能 是 因为 采 
用 的 超 晶 胞 不 够 大 或 者 原 胞 结构 的 原子 弛 豫 不 够 ， 计 算 中 要 避免 采用 有 虚 频 的 结果 ， DFPT 法 
得 到 的 Bi 色散 曲线 结果 并 没有 虚 频 ， 下 面 将 采用 DFPT 法 的 Bi 声 子 计算 结果 进行 热 中 子 散 射 数 
据 的 研究 。 


T x WK T L UW L i ; X WK T L UW L KUX 
(a) FP 法 (b) DFPT 方法 
4 FP 法 和 DFPT 法 计算 得 到 的 Al 声 子 色散 关系 


(a) FP 法 (b) DFPT 方法 
5 FP 法 和 DFPT 法 计算 得 到 的 Bi 声 子 色散 关系 


—s— WKress 
104 —e— Al Hawari | 
a ANTE 


00 
000 002 0004 0006 0.008 000 002 004 0016 


B kV 


图 6 不 同 模 型 Al 的 声 子 态 密度 的 差异 ”图 7 Bi 的 声 子 态 密度 计算 值 与 实验 值 的 差异 


Al 和 Bi 得 到 的 声 子 态 密度 计算 结果 与 文献 对 比如 图 6、 图 7 所 示 ， 由 于 ENDF8.0 和 FP 法 、 
DFPT 法 结果 的 声 子 态 密度 单位 不 统一 ， 本 文选 取 了 B 网 格 下 的 归 一 化 的 声 子 态 密度 进行 对 比 ， 
Al 原子 的 FP 法 和 DFPT 法 ， 以 及 ENDF8.0 的 声 子 态 密度 中 数据 基本 趋势 一 致 ， 从 对 数 图 的 结果 
可 以 看 出 ， 低 能 区 DFPT 和 ENDF8.0 的 符合 更 好 ， 相 较 于 Al 原子 ，Bi 的 声学 项 和 光学 项 分 隔 明 
显 , 计算 值 在 0.008-0.009eV 区 间 态 密度 值 为 零 ， 这 主要 归 因 于 Bi 的 结构 与 简单 的 立方 晶 格 的 差 
异 , 本 文 Bi 的 DFPT 计 算 结 果 相 较 于 Hawari 的 计算 在 声学 项 上 更 接近 W. Kress 的 实验 声 子 态 密度 
[由 ， 光 学 项 与 Hawari 的 结果 基本 一 致 ， 但 峰值 数据 偏 低 ， 并 且 有 更 小 的 展 宽 ， 两 个 计算 值 与 实 
验 值 均 稍微 有 偏差 , 出 现 偏差 的 原因 可 能 是 计算 中 原 胞 结构 的 原子 弛 豫 不 够 以 及 实验 本 身 测量 
的 误差 ， 整 体 来 看 ，DFPT 法 的 声 子 态 密度 产生 方法 更 灵活 和 可 靠 。 

3.3 热 中 子 散 射 截 面 数据 对 比 

分 别 基 于 FP 法 和 DFPT 法 得 到 的 Al 声 子 态 密度 , 利用 NJOY 完 成 了 Al 品 体 不 同 下 声 子 态 密度 
模型 热 中 子 散 射 截面 的 对 比 ， 如 图 8 所 示 ， 从 图 8(a) 可 以 看 出 ，FP 的 声 子 态 密度 得 到 的 AI 不 相 
干 非 弹性 散射 截面 比 DFPT 声 子 态 密度 结果 偏 高 ，DFPT 模 拟 结果 与 ENDF8.0 的 结果 更 接近 ， 主 
要 差异 体现 在 103 到 107 MeV 之 间 。 从 图 8 (b) 可 以 看 出 ， 相 较 于 FP 法 ， 采 用 DFPT 法 得 到 的 声 子 
态 密度 与 ENDF8.0 不 相干 非 弹性 散射 截面 和 相干 弹性 散射 截面 符合 较 好 ,， 并且 随 着 能 量 的 增加 ， 
FP 法 与 DFPT 法 的 截面 差异 越 大 ， 而 DFPT 法 和 ENDF8.0 基 本 一 致 ， 说 明 DFPT 法 的 声 子 态 密度 
与 ENDF8.0 采 用 的 声 子 态 密度 在 热 中 子 散射 特性 上 更 匹配 。 


1%- 
—s— DFPT 
i p pp 
4 _。 ENDF80 —e— ENDF8.0 
e per 
S g 
5 Él 
R B 
6 
a 
01 t : mr 10° r T T 
1E-10 1E-9 1E8 1E7 1E6 1E-5 10° 10’ 
Enerqy/MeV Energy/MeV 
(a) 不 相干 非 弹 散射 截面 (b) 相干 弹 散 射 截 面 


图 8 ”Al 不同 声 子 态 密 度 模 型 热 中 子 散 射 截 面 的 对 比 图 


hag 考虑 到 Bi 的 FP 声 子 色散 关系 有 虚 频 ， 本 文采 用 DFPT 法 得 到 的 Bi 声 子 态 密度 计算 了 Bi 原子 


Li 


300 有 下 热 中 子 散 射 截 面 ， 由 于 目前 国际 公开 核 数据 库 中 均 没 有 Bi 的 热 中 子 散 射 数据 ， 因 此 本 


文 将 计算 值 与 自由 气体 模型 和 实验 值 进行 了 比较 ， 如 图 9 所 示 ， 从 图 9 (a) 可 以 看 出 ，Bi 原 子 


由 气体 模型 在 低能 时 ， 特 别 是 能 量 在 布拉格 阔 值 以 下 ， 大 大 高 估 了 Bi 原子 的 中 子 总 截面 值 ， 这 


时 总 截面 主要 是 非 弹性 散射 截面 的 贡献 , 如 果 单 纯 使 用 自由 气体 模型 将 会 使 得 最 终 的 结果 与 实 


际 的 误差 比较 大 ， 热 中 子 过 滤器 Bi 对 热 中 子 的 过 滤 效 果 没 有 体现 出 来 。 从 图 9 (b) 可 以 看 出 ， 本 


文 的 计算 结果 与 实验 值 较为 一 致 ， 能 明显 体现 出 Bi 特定 阔 值 下 的 布拉格 效应 ， 同 时 也 能 看 出 ， 


在 布拉格 阔 值 以 下 ， 本 文 计算 的 非 弹性 散射 截面 与 实验 值 也 符合 较 好 。 因 此 ， 本 文 从 实际 需要 


出 发 ， 研 究 Bi 的 热 化 效应 是 有 必要 的 。 
0 0 e DJ.Hihes e 
四 到 -本 文 IT 作 è 
5 5 a 
£2 ~ 
H = 
ee I As 了 
GENTE Bic 
a spe 区 
w ype | 
E w e e 0 1w 105 104 103 10? 101 
人 EMNv 能 量 /eV 
(a) 与 自由 气体 模型 对 比 (b) 与 实验 值 进行 对 比 


Q 声 子 态 密 度 模 型 与 自由 气体 模型 及 文献 实验 值 的 对 比 图 


4 


结论 


本 文 以 A1H 和 Bi 金属 的 热 中 子 截面 制作 为 例 ， 分 别 基 于 冷冻 声 子 FP 法 和 密度 泛 函 微 扰 理论 


DFPT 法 , 采用 VASP 结 合 PHONONY 程 序 计 算 了 它们 的 声 子 色 散 关 系 和 声 子 态 密 度 , 基于 NJOY 


程序 
对 于 
于 Bi 


， 在 LEAPR 中 添加 Bi 相干 散射 处 理 ， 完 成 了 Al 和 Bi 的 热 中 子 散射 截 面 库 制 作 ， 结 果 表 明 ; 
Al， 相 较 于 FP 法 ，DFPT 法 的 声 子 态 密度 得 到 的 热 中 子 散 射 截面 与 ENDF8.0 符 合 更 好 ， 对 


， 采 用 DFPT 法 消除 了 FP 法 的 虚 频 现象 ， 得 到 的 热 中 子 散 射 结果 与 实验 值 符合 较 好 。 本 文 
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索 材 料 内 部 特征 机 理 的 角度 提出 了 一 种 更 基础 、 更 有 预见 性 的 制作 热 中 子 截面 的 方法 ,为 
新 型 反应 堆 的 核 材 料 热 化 机 理 莫 定 了 理论 基础 。 
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Application of First Principles in Calculation of Thermal 
Neutron Scattering Cross Section 


Lipeng Wang”, Xinyi Zhangl, Duoyu Jiang”, Tianliang Hul, Lu Cao! 


(1. Northwest Institute of Nuclear Technology, Xi’an, 710024;) 

Abstract: First principle is applied to the calculation of thermal neutron cross section of thermal 
neutron scattering materials in nuclear engineering field. In this paper, Al and Bi metals are used as 
examples to make thermal neutron cross section, respectively, based on the first principle, the frozen 
phonon method and density functional perturbation method are used. Using VASP and PHONONY, 
the phonon dispersion relation and phonon state density were calculated. Based on NJOY, Bi 
coherent scattering was added to LEAPR, and the thermal neutron scattering cross-section libraries 
of Al and Bi were generated. Results show that for Al, compared with the frozen phonon method, 
thermal neutron scattering cross section obtained by the density functional perturbation theory is in 
good agreement with ENDF8.0. For Bi, the imaginary frequency phenomenon of the frozen phonon 
method is eliminated by density functional perturbation method, and the results of thermal neutron 
scattering are in good agreement with the experimental results. In this paper, a more basic and 
predictable method for generating thermal neutron cross-section is proposed from the aspect of 
exploring the internal characteristic mechanism of the material, which lays a theoretical foundation 
for studying the thermal mechanism of the nuclear material of the new reactor. 

Keywords: First principle; Phonon density of state; Thermal neutron scattering cross section; 


Nuclear data generation. 
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